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Abstract 

We investigate the connection between the conservation of R-parity in supersymmetry and the Stueckel- 
berg mechanism for the mass generation of the B — L vector gauge boson. It is shown that with universal 
boundary conditions for soft terms of sfermions in each family at the high scale and with the Stueckelberg 
mechanism for generating mass for the B—L gauge boson present in the theory, electric charge conservation 
guarantees the conservation of R-parity in the minimal B — L extended supersymmetric standard model. 
We also discuss non-minimal extensions. This includes extensions where the gauge symmetries arise with 
an additional U(\)b-l &> U(l)x, where U(l)x is a hidden sector gauge group. In this case the presence 
of the additional U(l)x allows for a Z' gauge boson mass with B — L interactions to lie in the sub-TeV 
region overcoming the multi-TeV LEP constraints. The possible tests of the models at colliders and in dark 
matter experiments are analyzed including signals of a low mass Z' resonance and the production of spin 
zero bosons and their decays into two photons. In this model two types of dark matter candidates emerge 
which are Majorana and Dirac particles. Predictions are made for a possible simultaneous observation of 
new physics events in dark matter experiments and at the LHC. 
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I. INTRODUCTION 



R-parity is an important symmetry in supersymmetric theories (For a review see [1]). In 
supergravity theories Q2[| , over most of the parameter space of models consistent with the radiative 
breaking of the electroweak symmetry, the lightest neutralino is found to be the lightest supersym- 
metric particle, and this, along with R-parity (defined as R = (— i) 2S, + 3 ( B - L ) ? where S, B and L 
stand for the spin, baryon number and lepton number, respectively) and charge neutrality allows 
for the lightest neutralino to be a promising candidate for cold dark matter as suggested in [3|]. 

The question then, is, if indeed R-parity turns out to be a conserved symmetry of nature, how 
does such a symmetry come about, and how one may guarantee that it is conserved. It is known 
that the MSSM with the inclusion of a right handed neutrino, one for each generation, has an 
anomaly free U(1)b-l which can be gauged 1 . Of course, a U(1)b-l g au g e boson must grow 
mass otherwise it would produce an undesirable long range force. In the analysis that follows it 
is shown that a gauged B — L symmetry, where the gauge boson develops a mass through the 
Stueckelberg mechanism extending the Standard Model gauge group B4|] B5i|6|] preserves R-parity, 
i.e., R-parity does not undergo spontaneous breaking by renormalization groups effects under the 
assumption of universality of soft scalar masses, charge conservation and in the absence of a Fayet- 
Iliopoulos D-term. We will later refer to this model as the Minimal B — L Stueckelberg Extension 
of the MSSM. 

The fact that the minimal gauged B — L model proposed in this work preserves R-parity, with 
mass growth arising from the Stueckelberg mechanism, is in contrast to models with a gauged 
B — L where the symmetry is broken spontaneously and thus does not necessarily preserve the 
R-parity invariance. Thus the analyses of |M-Il2n show that R-parity symmetry, even if valid at the 
grand unification scale, could be broken by renormalization group effects 2 

We will first discuss the minimal (B — L) Stueckelberg extension of the Standard Model and 
of the minimal supersymmetric Standard Model (MSSM). In these extensions the Z 1 boson 3 is 



1 A gauged U (1)b-l arises naturally in GUT models such as SO(10) and Eg and in string models. 

2 For grand unified models where R-parity symmetry is automatic see II 1311 . For analyses where the spontaneous 
breaking of B — L occurs see Bl4ll5ll . for early work on t he sp ontaneous breaking of R-parity see lU6l4l9ll . For 
early analyses with R-parity and additional gauge fields see 120]. 

3 For recent dedicated work on heavy Z' B L physics see 112 1[ 12211 . 
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constrained to be rather heavy, i.e., it lies in the multi-TeV range and thus a direct detection may 
be difficult. However, this constraint is overcome in a U(1)b-l ® U(l)x Stueckelberg exten- 
sion, where U(l)x is the hidden sector gauge group. Here the Stueckelberg sector generates 
two extra massive vector neutral bosons, i.e., Z' and Z", one of which would be very narrow 
and could lie even in the sub-TeV region, and thus would be accessible at the LHC. The mod- 
els with massive mediators arise generally via mass mixing and kinetic mixing of Abelian gauge 



bosons [23 



■ 3311 . 041-13611 and the mixings are also the source of the so called dark forces [23, 25] 
- the mixings allow for a portal between the hidden (dark) sector via massive mediators R23l - l28ll 
(from which several components of dark matter can arise) and the visible sector where the states 
charged under the the Standard Model reside. Specifically, the class of models that we study here 
allows for two component (Majorana and Dirac) dark matter B7I1 . Such models with dark forces 
have received considerable attention in the context of the recent cosmic anomalies 071 . 139144111 ; 
for recent additional works on dark sectors see e.g. ll45l448ll. 

The organization of this paper is as follows : In sec. (fll]) we propose a U(1)b-l extension of the 
Standard Model via the Stueckelberg mechanism. In sec. (ITTl)) the B — L Stueckelberg extension of 
MSSM is introduced. In sec. (HVl) we outline the conditions for R-parity to be not spontaneously 
broken. In sec. (© we give a dedicated analysis of a U(l) b-l ® U(l)x extension of the MSSM 
via the Stueckelberg mechanism and show that the model naturally leads to a sharp Z' prime res- 
onance that can be seen at the LHC, and we analyze recent constraints from the Tevatron and the 
LHC. Here we also analyze the production and decay of new spin-0 particles. These scalars are 
the real parts of the chiral Stueckelberg superfields, where the imaginary part are the axions which 
are absorbed giving masses to the Z' and Z". In sec. (TVTl) we show that the model allows for two 
component dark matter, one consisting of neutral Dirac dark matter and the other of Majorana 
dark matter which produce a relic abundance consistent with WMAP [4911 . We also explore the 
detection possibility of dark matter with the recent limits set by the XENON and CDMS collabora- 

nn 

tions [50L I5111 which allows for direct detection constraints to be connected with the corresponding 
constraints on the Z' production at colliders. In sec. (IVIII) we give an overview as to how mod- 
els of spontaneous R-parity breaking can be distinguished from the R-parity preserving B — L 
extensions. Conclusions are given in sec. (| Villi) . 
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II. B — L STUECKELBERG EXTENSION OF THE STANDARD MODEL 



The B — L extension of the Standard Model provides a natural framework to understand the 
origin of neutrino masses since the three families of right-handed neutrinos, needed to cancel all 
anomalies, are used to generate neutrino masses. We first consider a U(1) B -l Stueckelberg exten- 
sion of the Standard Model with the gauge group 

SU(3) C ® SU(2) L <g> U{1) Y <g> U{1) B - L . (1) 

The mass growth for the U(l) b~l occurs via the Stueckelberg mechanism for which the extended 
Lagrangian is given by 

£ = £stT L + ^Yuk L + ^SM, (2) 
C^ L = - \c, u C^ - \{M BL C» + d,a)(M BL C» + Pa), (3) 
^Yuk L = Y v l L Hu R . (4) 

Here £sm is the Standard Model Lagrangian, l\ = (vl^l) and H = ia 2 H*. As usual, the 
Standard Model Higgs is H T = (H + , H°). The above Lagrangian is invariant under the B — L 
transformations 

6C li = d li X, 5a = -M BL \. (5) 
Added to the above is a gauge fixing term 

£gf = -^(W + M BL&) 2 , (6) 

so that the vector field becomes massive while the a field decouples. Additionally the interaction 
Lagrangian 

£st = 9blC^J BL) (7) 

couples the Stueckelberg field to the conserved B — L vector current J^ L . We note that the 
B — L gauge field has become massive with a mass M BL while maintaining the U(1) B -l 
invariance. Since B — L continues to be a symmetry even after the mass growth of the Z' its 
properties are rather different from the model where the B — L gauge symmetry is spontaneously 
broken through the Higgs mechanism. We will return to this in a later section. It is important 
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to mention that in this theory the neutrinos are Dirac fermions since there is no way to generate 
Majorana masses for right-handed neutrinos as in the canonical B — L model. This is a natural 
consequence coming from the Stueckelberg mechanism. 



In the above, a kinetic mixing term is possible leading to a generalized mass and 



anetic mix- 



520 where the 



ings for a massive U(l) which will then generally mix with the SM sector II23L 
hypercharge vector boson B mixes via both mass and kinetic mixings Il23n . One then diagonalizes 
the Stueckelberg mass and kinetic mixing together [^,| 53], [54], [55]. A further generalization to 
multiple U(l)s reads 

N v 1 N s N v 

" (8) 



2^2 



2 

/im ) i 



n=l 



m=l 



with N v Abelian vectors and N s axions, where B^ = and the other vector fields correspond 
to either hidden or visible gauge jwmmetries. Recent works with multiple additional U (l)s have 



indeed been discussed recently 13711521 Dq] . 140-14211 . Our analysis is restricted to non-anomalous 



extension of the Standard Model (for the anomalous case see e.g. D57H590)- In the analysis that 
follows we will assume the kinetic mixing is absent and instead investigate the pure Stueckelberg 
sector in the absence of mass mixing of the hypercharge B with the Stueckelberg sector. For recent 
works on the Stueckelberg Mechanism see e.g. Il43l445l. I60M63I1 and for early work in the context 
of strings see Q64£|. 



HI. B — L STUECKELBERG EXTENSION OF THE MSSM 

Here we construct the minimal U(l) b~l extension of the MSSM using the Stueckelberg Mech- 
anism. The supersymmetric extension of Eq. © is 

Cat = (M BL C + S st + S st )\ 2§ 2 , (9) 

where C = (C M , Ac, -De*) is the gauge vector multiplet for U(1)b-l, an d the Stueckelberg mul- 
tiplet is Sgt — (p + i&,ipst,Fs) where p is a scalar while a is the axionic pseudo-scalar. The 
supersymmetrized gauge transformations under the U(\)b-l are 

5blC = Cbl + Cbl , 5 B LS st = -MblCbl , (10) 
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where ( is an infinitesimal chiral superfield. Next we couple the chiral matter fields $j consisting 
of quarks, leptons and Higgs fields of MSSM. These couplings are given by 

^matter = * m e 2ffBi<3fliC $ m 1 9 2 § 2 (11) 

where Q BL = B — L and the sum is implicit over the chiral multiplets m and the interaction term 
of Eq. CO) couples the B — L vector field to fermions. We focus on the bosonic part of the extended 
Lagrangian which is given by 

^[0,1] = ~\c,uC^ - \Ml L Cl- l -{d,pf - ^Ml L p 2 

-\Dji\ 2 - g BL M BL p HQblU - \[Y,fl9BLQ BL f] 2 • (12) 

i 

The superpotential of the B — L extended theory is simply 

W = pH u H d + J2l Y uQH u u c + Y d QH d d c + Y e LH d e c + Y v LE u v c \. (13) 

gen 

Aside from the term Y v LH u u c Eq.(fT3l) is the superpotential of MSSM but without the terms that 
violate R-parity. 



IV. R-PARITY CONSERVATION 

As pointed out earlier, while the Stueckelberg mechanism gives mass to the B — L gauge boson, 
the Lagrangian of the theory, after the mass growth, still has a B — L symmetry and hence a 
conservation of R-parity (R = (— l) 25 + 3 ( B ^ L ) = (— \) 2S M. Here M denotes matter parity, which 
is +1 for Higgs and gauge superfields, and —1 for all matter chiral superfields). This conservation 
of R-parity in the minimal B — L Stueckelberg extensions is in contrast to models where the B — L 
gauge symmetry is broken by a Higgs mechanism and where in general the mass growth of the 
B — L gauge boson could break the B — L symmetry and thus R-parity invariance is also lost. For 
example, for the model of Eq. (fT3~l) . a VEV growth for the scalar field in the uf multiplet will break 
B — L invariance and generate a mass for the B — L gauge boson. However, a VEV growth for 
u° also violates R-parity invariance which then removes the neutralino as a possible candidate for 
dark matter. Specifically, for example, in Eq. (fT3l the VEV growth of vf generates the term LH U 
which breaks R-parity. However, in the minimal B — L Stueckelberg extension of MSSM even 
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after the mass growth of the B — L gauge boson R-parity is maintained and the R-parity violating 
interactions such as LH U , LLe c , QLd c , u c d c d c are all forbidden in the superpotential. 

A. Scalar Potential and R-Parity Conservation 

We wish to show here that with a Stueckelberg mechanism for mass generation the B — L 
symmetry not only remains unbroken at the tree level but further that this invariance is not violated 
by radiative breaking in the minimal model. We give now the deduction of this result which is 
rather straightforward. We exhibit below the potential including just one generation of leptonic 
scalar fields in the model consisting of p, z>, e, e c , v c (An extension to 3 generations is trivial). 
Assuming charge conservation so that (e) = = (e c ), etc., and including soft breaking, the 
potential that involves p, v and v c fields is 

Vst-BL = l -{M 2 BL + m 2 f )p 2 + M 2 u^ + M 2 c ^u c 

+ ^ {v c] v c - + 9bl M BL p (u c ^ c - tfv) , 

+ \Y u \ 2 {\Hy\ 2 + \uH Q u \ 2 + \uu c \ 2 ) 

- \(9 2 + 9' 2 )(\H° U \ 2 - \H°\ 2 )\v\ 2 + \(g 2 + g> 2 )\v\* 

+ (-p*Y u H° d *uu c + A v Y v VV c Hl + h.c.) (14) 

where we have used Qbl(z) = Qbl{v) = —1 and where m p , Mp, and M^c axe soft masses. The 
relevant part of the potential is then 

V = } Vst-BL + Vmssm, (15) 
•'—'gen 

and where as is familiar 

Vmssm = W + m 2 H J\H° u \ 2 + (\p\ 2 + m 2 Hd )\H° d \ 2 - (BpH° u H° d + h.c.) 

+l(9 2 + 9 ,2 )(\H° u \ 2 -\H° d \ 2 ) 2 - 06) 

We begin with universal boundary conditions for the RGEs. We note that the RG evolution for M g 
and Mjj are identical since SU (2) L ®U(l) Y symmetry is unbroken down to electroweak scale. If 
Mf turned tachyonic it would lead to VEV formation for the field e violating charge conservation 
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and thus we disallow this possibility. Since v and e lie in the same SU (2) l multiplet the same 
holds for the v field, i.e., it too does not develop a VEV. This can be seen from the one loop RG 
sum rule connecting the sneutrino v mass and the selectron mass 

Ml - Ml = cos(2/3)M^ + <% e) (17) 

where b 2 v e is difference of the mass squares of the fermions (and is essentially negligible compared 
to W mass term the largest of which occurs for e — > t which is still negligible). Thus the right 
hand side of Eq.fPTl) is positive definite for any range of tan/3 in the perturbative domain in the 
RG analysis. As a consequence, if the mass square of e does not turn tachyonic, this also holds for 
the mass square of v and (D) = 0. Thus with (e) = = (p) = (e c ) , and integrating on the p field, 
we get the following potential for v c 

V " = M ^ + 2(^t% { ^ C)2 + ' y " |2|Wy|2 ' <18) 
The last term above is negligible in size compared to the other terms since it involves the Yukawa 
Y v . Thus the coupling between this sector and the MSSM sector via the H° field is negligible. 
Now in the RG analysis there are no beta functions to turn M~ c negative and the quartic term is 
positive definite so the potential is bounded from below. Consequently the potential cannot support 
spontaneous breaking to generate a VEV for the field v c and thus (i> c ) = 0. Further, the extrema 
equation for p gives 

<'> = -]|^<^-^> = °' ™ 

and since (z>) = = (z> c ), one also has (p) = 0. Thus there is no spontaneous symmetry breaking 
in the system and the B — L and consequently an R-parity is preserved. We add that the situation 
here is rather different from the Stueckelberg extensions introduced in [|4j-l6|] where p receives a 
non- vanishing VEV. In |4-6fl, a non- vanishing VEV for p would arise due the Stueckelberg sector 
mixing with the U(l)y sector of MSSM. In contrast in the minimal B — L extension analyzed 
here there is no mixing with the U(l)y sector, and thus there is no VEV growth for p. Thus the 
entire mass growth in the £ 7(1 ) b—l sector occurs via the Stueckelberg mechanism. If we include 
a Fayet-Iliopoulos D term [65I] then effectively the potential for v c is replaced by 

V - = M ^ + 2(Ai L '+W ^ + « 2 + WW < 20 > 



For the case when £ is negative a VEV growth for v c is possible and R-parity can be broken 
spontaneously. While an FI D-term naturally arises when the U{\) is anomalous the inclusion of 
an FI term for a non-anomalous U (1), which is the case we discuss, is superfluous, and we exclude 
it from the minimal model. Therefore, it is apparent that R-parity is always conserved within the 
minimal Stueckelberg B — L extension of the MSSM. 

The analysis above follows with (minimal) universal boundary conditions on the soft scalar 
masses. However, since the nature of physics at the Planck scale is still largely unknown one 
should consider non-universalities as well. In this case one will have additional contribution to 

— ii — i r~i 

the mass squares of scalar masses [|66i|67|]. The analysis of [68] considers a contribution to M~ c 
arising from Tr{Q BL m 2 ) with 

Sbl = Tr(Q BL m 2 ) = 2(M§ - Mf) + (M| c - M%) + (M 2 C - M| e ), (21) 

under the constraint Tr(Ym 2 ) = 0, where 

S Y = Tr(Ym 2 ) = M 2 H2 - M 2 Hi + ^(M| - 2M 2 C + M% - M\ + Mf c ). (22) 

gen 

With the universal boundary conditions for only each family one has Sbl = 0. This can be 
achieved in minimal supergravity models where all scalars have the same soft mass term, or in 
50(10) or E 6 scenarios where the boundary conditions tell us that all sfermions of one family 
should have the same soft mass term. However, with non-universal boundary conditions one will 
have in general Sbl ^ 0. With inclusion of Sbl one could in principle turn M 2 C negative. Such a 
situation is achieved with inclusion of specific constraints in the analysis of [68]. However, such 
constraints are not generic and the positivity M~ c may still be broadly valid even with inclusion of 
non-universalities of soft parameters. 



Now there are stringent bounds on an extra B — L type gauge boson. One finds [69] 



M z '/gBL > Mbl ~ 6 TeV, (23) 

which implies that for gBL ~ 1 the B — L type Z' boson lies in the several TeV region. With a Z' 
of this mass scale, detection at LHC-7 may be difficult, both because of energy considerations and 
luminosity. Further, with the constraint as given by Eq.(T23l some of the other phenomenological 
implications of the model associated with the spin and spin | sectors will also be difficult to 
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test. In what follows, we uncover a model which maintains the strict R-parity invariance of the 
minimal Stueckelberg B — L extensions, even after mass growth of the B — L gauge bosons, but 
with testable implications that are far more rich. 



V. 17(1) B—L <S> U(l) x STUECKELBERG MODEL 



As indicated in the last section, the Z' boson of the minimal B — L model may be difficult 
to detect because of its heavy mass. We consider now an extension of the model of the previous 
section which overcomes this constraint and produces a Z' which is much lighter but still has 
B — L interactions with matter. This extension includes a hidden sector U(l)x which is anomaly 
free but allows for a mixing between the visible and the hidden sectors. The extended gauge group 
reads: 

SU(3) C ® SU(2) L ® U(1) Y ® U(1) B -l <8> U(l) x ■ (24) 
Thus we have Stueckelberg mass growth in the Abelian sector via the interaction 



-St 



d 2 9d 2 9 [(MiX + M' 2 C + s + sf 
+(M[X + M 2 C + S' + S') 2 } , 



(25) 



where the model is invariant under the extended gauge transformations 

8 X (X,S,S',C) = (e x + ex,-M ie x,-Miex,0) 
5 BL (C,S,S',X) = (e BL + e BL ,-M^ BL ,-M 2 e BL ,0) 



(26) 



where ex,BL are infinitesimal chiral superfields. One can compute the mass matrix for the U (l)x 
and the U(1)b-l gauge vector bosons by going to the unitary gauge which in the basis X^, 
gives 



^[spin 1] 



M\ + M' 2 Mi M' 2 + M[ M 2 
/', + M[M 2 M 2 + M' 2 



(27) 



Here M[, M' 2 are the mixing parameters and in the limit that M[, M' 2 — y we have that the masses 
of the Xp, C M bosons are M 1: M 2 . The diagonalization gives us two massive vector bosons which 
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// T(Z> -> ff)/a BL M z , T{Z" -> ff)/a BL M z „ 



sin 2 9 BL /3 cos 2 6bl/3 

vt-vt. sin 2 6^/3 cos 2 9 B l/3 

qq(q + t) f s sin 2 6 BL /9 f s cos 2 9 BL /9 

tt f s f tiZ , sin 2 BL /9 f a f t , Z n cos 2 BL /9 

TABLE I: The decay widths of the Z' and of the Z" bosons into leptons and into quarks in the U(l)x <8> 
U(1) B -l Stueckelberg model where a B L = g BL /A7r and f s = (1 + ^) and for 1/ = Z', Z", one has 



Av = (i + 2^)(i-^). 



we may call Z', Z" where 

= cos 8 slZ'^ + sin 8 slZ'^, 
Cf, = -sn\d BL Z^ + cos 8 BL Z'^. (28) 

We consider now the case of small mixing, i.e., M[, M' 2 <C Mi, M 2 which implies tan^i <C 1. 
For small mixings the Z' boson lies mostly in the hidden sector with a small component propor- 
tional to tan 8 B l in the B — L sector while the opposite holds for Z" . Here Z" lies mostly in the 
B — L sector with a small component proportional to tan 8 B l in the hidden sector. 

Since lies in the hidden sector and has no couplings to the visible sector matter, the only 
couplings of Z' ', Z" to the visible sector arises because of the couplings of C M to the visible sector 
matter. Using the couplings of C M one finds the couplings of Z' and Z" to the fermions (/j) to be 
of the form 

C Z i,z» = ij il 11 9blQ Bhf i) [— sin 8 B lZ'^ + cos8 BL Z;}. (29) 
In the context of Eq.d29l) the constraint of Eq.d23l) gives two separate conditions, i.e., 

M z >/g BL > sin8 BL x (6 TeV), 

M z n /g BL > cos 8 BL x (6 TeV) . (30) 

It is clear that the constraint on the Z' is now considerably weakened relative to the constraint of 
Eq.d23l) if the mixing angle 8 BL is small and one can have 

M z , < 1 TeV, Stueckelberg C/(l) B -x ® U(l) x . (31) 
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FIG. 1: Upper panel: An exhibition of a(pp — > Z') • Br(Z' — > e + e~) vs the mass of the Z 1 resonance in 
the Stueckelberg U(1)b~l &> U(l)x extension of MSSM at the Tevatron. Here qbl = 0.35 and sin#BL 
takes on the values 0.01, 0.05 from the bottom to the top curves in the plot. The analysis assumes that the Z' 
decay into the hidden sector is suppressed. Lower panel: The same analysis at LHC-7 with sin 6bl taking 
on the values (0.01,0.02,0.03,0.05) from the bottom curve to the top in that order. 

However, Z" is still heavy since cos 9 bl ~ 1 for small 6 B l- In Table® we give the decay 
widths of the Z' and Z" bosons into leptons and into quarks. The relative strength of the Z' 
decay into quarks and leptons provides a distinctive signal for this model. Thus, for example, the 
ratio of the branching ratios of Z' into charged leptons vs into quarks (except into tt) is given by 
BR(Z' £ + £-)/BR(Z' -> qq) = 6/(5(1 + ajir)). Further, in this model the decay width of 
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Stueckelberg u(i) BL x u(i) x 
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FIG. 2: Exhibition of a 500 GeV Z' resonance in the Stueckelberg U(1) B -L ® U(l) x model at LHC-7 
with a variable luminosity from 5fb _1 to lOfb^ 1 with a Pt cut on leptons of Pt > 30 GeV. Currently the 
LHC has analyzed ~ lfb -1 of luminosity. For a Z' resonance of 500 GeV with 8bl = 0.05 and qbl ~ 9y 
the LHC would need about 5fb _1 to begin to see any Z' effect. With a very optimistic 20fb _1 , the Z' signal 
will be strong and Z' should be visible with the mixings and masses of the size discussed. 



the Z' and Z" are related by 



(32) 



r(z"^j: i fifi) ™ 

Eq.(l32l) implies that for the Z' mass in the sub TeV range, and the Z" mass in the range above 
6 TeV, and tan 6 BL 1 consistent with Eq.(l23l). the ratio of the decay widths of Z' vs of Z" can 
be vastly different, i.e., a decay width of Z' in the MeV range vs the decay width of Z" in the 
hundreds of GeV range. Thus while the Z' will be a very narrow resonance, the Z" will be a very 
broad resonance. 

It is also instructive to check the contribution of the new interactions to the muon anomalous 

O 

moment which is now measured very accurately [70] so that the current error in the determination 
is given by A(g M — 2) = 1.2 x 10~ 9 . The contribution of the Z' and of the Z" bosons to the 
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anomalous moment is given by 



24vr 2 



sin 2 6 BL cos 2 6 B l 



Ml„ 



(33) 



Using the LEP constraint of Eq. (|30l) one finds that the contributions of the new interactions is 



A(o, - 2) < 



m. 



127T 2 M 2 L 



(34) 



and a substitution of M B l ~ 6 TeV gives a rather small contribution, i.e., A(g /t — 2) < 0(1) x 
1CT 12 . Remarkably in this case the LEP constraint of Eq.(l30l) is stronger than the constraint arising 
from the very precise measurement of g M — 2. 



A. Production of Vector Resonances 

The fact that the Z' boson could have a low mass has important phenomenological implications. 
From Table® we note that the decay width of the Z' boson is proportional to sin 2 9bl and since 
sm9 BL is small, the decay width is relatively small, i.e., with the mass of the Z' in the sub TeV 
region, its decay width would be in the MeV range and thus the Stueckelberg Z' is a very narrow 
resonance. A narrow resonance of this type should be testable in collider experiments much like 



the hypercharge Stueckelberg Z' on which the D0 currently has experimental bounds [71]. Further, 
the decay of the Stueckelberg Z' into leptonic channels will be much more than in the hadronic 
channels because the branching ratios are proportional to (B — L) 2 . Thus one can discriminate 
a B — L Stueckelberg Z' boson by a study of its branching ratios. Such a resonance could be 
produced in the Drell-Yan process at the LHC and the Tevatron via 

VVipP) — ► Z' — ► QQ- (35) 

In Fig. (OQ) we show the predictions for the Z' cross section times the branching ratio into e + e~ 
in the U(\)b-l ® U(l)x extension of the Standard Model. Cross sections and event rates are 



calculated by implementing the couplings into PYTHIA and PGS [|72i 173(1 • The bottom panel 
shows the limits on the production cross section for a(pp — > Z' — > ee) at y/s = 7 TeV with the 
recently released ~ 1 fb _1 run [ 74I1 . For these curves we take g B L to be have the same value as 



the hypercharge gauge coupling g Y and we let sin6> BL run from 0.01 to 0.05. The cross section 
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for other values of the product g B zsm9 B L can be estimated by the scaling in the cross section 
which at the Z' resonance scales like g 2 BL sin 2 9 B l- The top panel give s a similar analysis for 
the Tevatron using the D0 data with 5.4/fb of integrated luminosity 117 IN . From the analysis of 
Fig. £[]) we observe that at present the Tevatron bound is about as strong as the present LHC bound. 
However, the LHC will surpass the Tevatron very soon. Indeed, the Z' produced in the model can 



exist with a much lower mass 



75 



76|] than the Z' models presently excluded by ATLAS [LZJ 



and CMS flVVh . In Fig.© we display the number of events as a function of the di-lepton invariant 
mass. Here one finds that with an optimistic choice of an integrated luminosity of 20fb~ 1 the 
number of dileptonic events in excess of 30 in the peak mass bin and should be visible. Thus a 
Z 1 mass of 500 GeV with a mixing angle 9 BL = 0.05 and Qbl = 9y is a promising candidate for 
discovery. 



B. Production and Decay of the Scalars p and p' 

In addition to the Z' phenomenology there are other sectors where new phenomena can arise. 
One of these relates to the scalar components p x and p B L of S + S and of S' + S' that remain in 
the bosonic sector after Z' and Z" gain mass by the Stueckelberg mechanism. These fields mix 
with the D-terms so that one has the following set in the Lagrangian 

P bl{M x D bl + M' 2 D X ) + p x (M[D BL + M 2 D X ). (36) 

Elimination of the D-terms gives the following mass matrix in the p x and p BL basis 

\m 2 + M' 2 + ml M : .\l{ + M 2 M' 
M l M[ + M 2 M! 2 M[ 2 + M 2 + m 2 BL 

where we have also included the soft contributions to masses for p x and p B L- We note that the 
structure of the spin zero mass squared matrix given by Eq.(l37l) is different compared to the mass 2 
matrix given by Eq.(l27T). The reason for this is that while the vector mass squared matrix arises 
directly from the Stueckelberg term Eq.(|25l), the mass squared matrix of Eq.(f37]) arises from the 
mixing given by Eq.(T36l). The mass matrix of Eq.([37]) gives two mass eigenstates p and p' with 
eigenvalues M p and M p i. The mass parameters M[, M' 2 can define the mixing and when the mixing 
is small, M 2 M\ + m\ and M 2 , ->■ Mf + m 2 BL . With M p in the sub TeV range M p > may have 
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a mass in the several TeV range. These mass eigenstates are admixtures of px and pbl so that 
px = cos 0' BL p + sin 9' BL p' and p B L = — sin 8' BL p + cos 9' BL p' . For the case when the soft terms 
are absent, the eigenvalues of the mass squared matrix of Eq.dTTT) is are identical despite the very 
different looks of the matrices of Eq.(|2Tb and Eq.dTTl). This can be seen by the following unitary 
transformation 

[/ t M [ 2 spinl] [/ = M [ 2 spin0] , (38) 
where the unitary matrix that connects the spin 1 and spin matrix is given by 

/cos 5 sin A nf= ff-ff , m 
\-smC coaCj Mi + M 2 

This result shows that the eigenvalues for the matrices M? in ^ and M? ia Q , are the same in the 
limit of vanishing soft masses for the scalars. Now it is assumed that all the matter fields in 
the visible sector do not carry any U(l)x quantum numbers, i.e., Q x = for quarks, leptons 
and the Higgs fields. Further, following the analysis of Sec. (4), it is straightforward to establish 
that the quartic term (z> c ^z> c ) 2 has a positive co-efficient in the scalar potential. Thus once again 
since there are no couplings in the model to turn M? c negative, there is no spontaneous violation 
of R-parity also in this extended model while the B — L gauge boson develops a mass via the 
Stueckelberg mechanism. 

From the discussion preceding Eq. (|38l) , it is clear that the field px has no coupling with the 
visible sector while p BL has couplings of the form g BL Mp BL fiQ BL fi. One then has the following 
interactions of p and p' with sfermions 

C pN = - smO , BL g BL Mj}Q BL f i p + cos6 BL g BL Mj}Q B J iP '. (40) 

Eq.(l40l) allows the decay of the p(p') via its couplings to the sfermions. If kinematically allowed 
p(p') will decay into leptons + E™ tss or into jets + E™ tss where E™ tss contains at least two neu- 
tralinos x° (here x° is the lightest neutralino (LSP) of the U (1) B -l ® U (l)x combined sector - see 
sec. (fVTl) ). However, an interesting situation arises when the mass of p(p') is smaller than 2M x o . In 
this case p(p') cannot decay into the final states with 2 LSPs and only the decays into the Standard 
Model particles are allowed. Such decays can occur via loops and the final states will consist of 
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FIG. 3: Diagrams giving rise to the production of the Stueckelberg scalars, p, p' at the lowest order. 



gg, fifi, WW, ZZ, 7Z, 77. There are many diagrams that contribute. The dominant one relevant 
to the model we study here with real scalars p and p' are the gluon fusion diagrams (see Fig. (3)). 

From Eq.(|40l) the interactions of p and p' to the mass diagonal squarks are given by the follow- 
ing interaction 



£( P ,p'W<i = -# P Micos(2%) (qliQup - ql&ip) - g P M 1 sm(26q i ) (ql&iP + qhqup] 

+ (p^p\- sin 9 BL ^ cos e BL ). (41) 



with the B — L dependance encoded via 



9 P = -^9bl sin 0' BL , 



(42) 



and where i runs over the squark flavors. Now while the p, p' vertices allow couplings with squark 
mass eigenstates, where the two states couple to are either the same state or different states, the 
gluino only couples to squark states, where both states have the same mass. Thus in Eq. (l4T|) only 
the interaction terms proportional to cos 20^ enter in the gluon fusion diagram. As such, the decay 
width of the p to gluons is given by 

2 



T(p gg ) = ^ 3 \ 1 

KH yyj 5127T 3 



E (-l) 1+a cos(2% 



Li(r ai ) 



a=l,2:i 



2 

m a 



(43) 



with r ai = M 2 / (Am^ ai ) , and Li(r) is a loop function defined by [|78|l 



Li{r) = r 2 [r - f(r)] , f(r) 



arcsin 2 (\/r) 

_I (\ oz l+v 7 !^ 7 



— m 



r < 1 
r > 1 



(44) 
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As a consequence of the symmetry of gauge interactions one also has 



M 3 , 

rV ^99) = ^ cot 2 ff BL Y{p gg). (45) 



Further the partonic production cross section of p is given by 



256vrM p 4 



^ (-l) 1+o cos(20 S( )r ai L 1 (r o 

a=l,2;i 



5(l-M 2 Js). (46) 



The hadronic production cross section relevant to the search for p at the LHC is a(pp — > p) and is 
given by a convolution with the parton distribution functions for the gluon, which at leading order 
in the narrow width approximation is given by 

dL pp 

a{pp -> p)(s) = r p —^a{gg -> p). (47) 
Here is the center-of-mass energy, r p = M 2 / s, and dJJ^Jdr is given by 

/ —f g /p{x,Q)f 9 /p{~,Q), (48) 
j T x x 



dr 

where / s / p is the parton distribution function for finding the gluon inside a proton with momentum 
fraction x at a factorization scale Q. A numerical analysis shows that a(pp — > p) can lie in the 
range 0(1000) fb in the most optimal part of the parameter space for producing the p. 

The final decay modes of the p can produce visible signatures at the LHC, and branching ratios 
will generally be different from the Standard Model Higgs hsM- Thus hsM has both tree level 
decays into the final states bb, rf, cc as well as decays via loop diagrams into gg, WW, ZZ, Zj, 77. 
For a Higgs boson mass of 100 GeV, dominant decays modes are the tree level decay modes with 
bb decay being almost 80%. Among the loop decays the dominant decay is gg and sub-dominant 
decays are WW (off shell) and 77 at a Higgs mass of 100 GeV. Now suppose the tree decays 
of the Higgs were suppressed, then the decay of the Higgs to 77 will have a branching ratio of 
~ 2.5 x 10 -2 . The decay of the p parallels this case since there are no tree decays of the p. In the 
analysis below we will use the above branching ratio to get an approximate estimate of 77 event 
for the p decay. An analysis of pp — > p at the LHC at y/s = 14 TeV is given in Fig.©. One finds 
that the cross section at M p = 100 GeV for the maximal case with (0' BL , g BL ) = (n/2, 1.2) is 
~ lOOfb. At 200 ftr 1 of integrated luminosity at the LHC at ^/s = 14 TeV, one will have 2 x 10 5 
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FIG. 4: A display of the production cross section a(pp — > p) from gluon fusion as a function of the mass 
of the Stueckelberg scalar p at the LHC with yfs = 14 TeV for several combinations of 6' BL and gsi for 
the case which maximizes the production for the MSSM sector, | cos 26*^ | = 1, i.e. \A t \ tan = \p\ . From 
bottom to top the curves have {0' BL ,gBL) = (vr/6, 0.65), (-7r/2,0.65), (vr/2, 1.2), where the top curve is 
close to the theoretical upper limit on the production. The kink appears at the point where M p j ftm^ ) = 1 
and the analysis has the other squarks and the gluino much heavier than the lighter stop. 

p events when M p = 100 GeV. Using BR(p — > 77) = 2.5 x 10~ 2 one finds ~ 5000 77 events 
before kinematic and efficiency cuts. We note that the photons coming from the 77 signal will 
be monochromatic carrying roughly half the mass of the decaying particle. Thus the 77 signal 
arising from the decay of the p would be distinguishable from the 77 signal from the Higgs decay 
if the masses of the two are significantly separated. A p mass of 100 GeV would imply a Z' mass 
of also 100 GeV assuming no soft terms in the p sector. A Z' mass of 100 GeV is consistent 
with the current data if either the mixing angle 9bl is small or the Z' decays dominantly into the 
hidden sector (see Sec. dVIBI) ). We note also that while the mass of the Z' and the mass of p are 
the same in the absence of soft breaking terms for p, the couplings of the Z' to fermions and of 
p to squarks can be of very different sizes. This is apparent from Eq.(l39l). Hence the possibility 
arises of being able to discover both the p and the Z' . However it is also quite possible that only 
one resonance may be visible depending on the overall size of the Stueckelberg masses and the 
individual couplings of the two states. 

20 



LHC Production : Stueckelberg Scalar 




The production cross section for pp — > p, p' bears resemblance to the analysis of Ill2h and is 
closely related to canonical Higgs production (see e.g. [78, 3]) but is restricted by the form of 
the couplings as given in the B — L Stueckelberg extension. We add that recently several models 
with scalars have been studied in the literature which can produce large production enhancements 
relative to the SM higgs production (see e.g [80-83]). The production of p does not receive 
enhancements of the size studied above, but nevertheless does produce event rates that can be 
measured at the LHC-14 with larger luminosity as was detailed above. 

We note that very recently the LHC has put new constraints on the allowed mass of the Stan- 
dard Model Higgs Boson h SM . Preliminary analyses based on those reported at EPS 2011 and 
at Lepton-Photon 2011 [8^] imply that the SM Higgs boson has a mass below ~ 145 GeV. The 
above result is compatible with the SUGRA models which typically indicate a Higgs mass below 
~ 140 GeV. Because the production of p relative to the hsM differs markedly via their cou- 
plings, as discussed above, the production of the two fields could be distinguished with sufficient 
luminosity. This is possible if the h$M resonance and the p resonance are sufficiently separated 
in mass. In addition, because the production of p is weaker than hsM, the golden channels such 
as ZZ, WW remain available where hsM has been ruled out to have such a mass. Searches for 
M p ~ (200 — 500) GeV will however have to wait for upgraded luminosity at the LHC. 



VI. NEUTRAL DIRAC AND MAJORANA COMPONENTS OF DARK MATTER 



A. Majorana Dark Matter 

The U(\)b-l ®U(\)x Stueckelberg extension of MSSM have new implications for the nature 
of dark matter. Specifically in the neutralino sector we have in addition to the MSSM neutralinos, 
extra gauginos and stinos, where the stino is the analogue of the higgsino. Thus from the gauge 
supermultiplets X = (X^, Ax, Dx) and C = (C M , Ac, -De) we can construct two gaugino states 
which we label as Ax,Abl- Similarly from the chiral multiplets S + S and S' + S' we can 
construct two higgsino states ipg, ijjg>. These four neutralino states in the Stueckelberg sector have 
no mixings with the MSSM neutralinos. Thus the neutralino mass matrix in the U(l)x ®U(1) b-l 
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extension of MSSM has the form 



neutralino 



Specifically the neutralino mass terms in the U(l)x <8> U{\)b-l sector are given by 



M st 


04x4 


04x4 


A^MSSM 



(49) 



£ mass = _ Z T Mgt Z (5Q) 

Z T = (i>s, 4>s>, Ab-l, Ax), (51) 

where the 4 x 4 sub-block of the U(1)b-l <8> U(1)x sector has the form (omitting for simplicity 
the soft terms) 

I 2x2 m \ ( M ± M>\ 

M st =\ , m = \ . (52) 

I m T 2x2 / \M[M 2 

\ / 4x4 \ / 2x 2 

We can diagonalize the neutralino mass matrix in the U(l)x ®U{\)b-l sector by an orthogonal 
transformation Z = OX so that 

X T T Ai st OX = diag(m x o,m x o,m x o,m x o). (53) 

Now the generalization of the matter Lagrangian reads 

^matter = $ m e 2 ^ L Q BL C + 2 gxQx X ^ ^ ^ (54) 

and gives a coupling of the type A B _ L f L fl, where are a chiral fermion and a chiral 

scalar, which leads to couplings of the Stueckelberg sector neutralinos with matter of the form 
(k = 5 — 8). Thus we note that even though the neutralino mass matrix does not have a 
mixing between the MSSM and the Stueckelberg sectors, the neutralino in the Stueckelberg sector 
can decay into the least massive supersymmetric particle (LSP) which may lie in the MSSM 
sector. The way it occurs is as follows: The neutralinos xl (k = 5 — 8) have fermion- sfermion 
interactions as indicated above, while the neutralinos in the MSSM also have similar type 
interactions. If the mass m x o > m x o we will have decays of the type xl ~^ fifi ~> fifiXi + ' ' ' • 
Thus there is only one stable Majorana supersymmetric particle in the combined MSSM and 
Stueckelberg system. On the other hand if, for example, xl is the lightest neutralino then the LSP 
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will lie in the Stueckelberg sector. In this case the Xst = X5 would be a dark matter candidate 
(the notation st denotes stueckelberg and does not imply preference of the stino component over 
the gaugino component). Its properties are expected to be similar to those of the bino LSP of 
the MSSM. For the case of a thermal relic, the annihilation of Xst will occur via the t-channel 
squark exchange so that (dropping the superscript from here on) Xst + Xst fifu as wells 
Xst + /mssm — > SM SM', x s t + Xmssm — > SM SM' where the last two cases indicate that the 
the coannihilations will generally occur I125n . D9IU 401 (for a review see [85]). For the direct 
annihilations, unlike the annihilation of MSSM neutralinos, there are no direct channel Z or 
Higgs pole exchange diagrams and consequently final states such as WW, ZZ, ZH, HH are 
absent at the tree level. For the case of co-annihilations this is modified. If p is of low mass, as 
discussed in the previous section, the stop should be relatively light to accommodate a signal of p. 
In this case the relic density can be satisfied via stop co-annihilations [96]. 



Next, we discuss the the direct detection of Xst- Specifically there are no t-channel Higgs 
or Z pole exchange contributions to the direct detection rates for this case at the tree level. As 



pointed out in Ref. [|86L 



870 h is important to include contributions arising in the spin independent 



scattering cross section from the twist-2 operators 



f p /m p 3 Yl fJr q + £ -(q(2) + q(2))g^ + ... 



(55) 



q=u,d,s 



q=u,d,s,c,t 



where the additional terms are suppressed and q{2), q{2) are matrix elements and are given in [87] 
Specifically g!p is given by, in the limit of massless quarks 



4" 



M. 



Xst 



a 2 q + b 2 q 



M 2 ) 2 



(56) 



where a 2 + b 2 q = q % 1 Q% L = g 2 PL /9. In addition, there are terms of size J2„ =u d s f q fr q (where 



f q , f Tq are given in |86j,|88L|89J]). Here terms in f q that are proportional to a 2 + b 2 are suppressed 



by a factor of 4 relative to g q ^ [87]. Terms in f q also contain a 2 — b 2 oc g 2 BL Q 2 BL sin 29 q and 
are ultra suppressed by the smallness of the squark mixing angle. For the case when the M Xst is 
relatively close in mass to m,, up to correction in the light quark masses, there is an enhancement 
in the SI cross section ll87ll . Utilizing this effect, for mass splitting of order 30-100 GeV, one easily 
sees detectable size SI cross sections for squark masses that are in accord with LHC limits (see 
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Stueckelberg Majorana Dark Matter 
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FIG. 5: Spin independent Xst neutralino-proton cross section vs the Stueckelberg neutralino mass for the 
case when the Stueckelberg neutralino is the LSP. Exhibited is the spin independent cross section for several 
combinations of A ~ rriq — M Xst (in units of GeV). The current limits from XENON- 100 are also exhibited. 



Fig.©). At even smaller mass splittings, the models are constrained by XENON. We have verified 
using micromegas ll9oll that the small mass splitting between the LSP and the squarks can lead 
to cross sections of the size we find. In this case the relic density can be brought in accord with 
WMAP from co-annihilatons. In particular the squarks in the initial state annihilations play a large 
role in reducing the relic abundance. There is also mixing that derives from rotating between the 
chiral fermion in the Stueckelberg multiplet. We consider the optimal case where in the mass 
diagonal basis, the lighter of the two mass eigenstates is the one which couples via the larger 
mixing. Thus we have taken the mixing in the gaugino stino sector cos 9 Xst — > 1, and have fixed 
9bl = 0.65 in Fig.©. The result of a large scattering cross section does require an LSP above 



around (500-600) GeV to be consistent with the current limits from the LHC [|91 



95] 



24 



B. Dirac Dark Matter 



Additional matter fields in the form of Dirac fermions (and their super symmetric counter parts, 
two chiral scalars) can exist in the U (l)x sector which have only vectorial couplings to the gauge 



field X M and a mass for the Dirac fermions can be generated via terms in the superpotential [37]. 
As seen already, after mixing of the B — L gauge field C M with the field X M , two mass eigenstates 
Z' and Z" arise in the mass diagonal basis each of which have B — L type couplings with the SM 
fields. In addition, the interaction of the dark sector Dirac field with the Z\ Z" is given by 

C D = D-fiCvoZ'p + C Z n D Zl)D. (57) 

The interaction vertices with the Dirac particle (D) with the visible sector quarks and leptons enter 
through the vector mixings so that 

C Z 'D = gxQxcos9 B L, C Z "d = gxQx sin B l ■ (58) 

The dark sector Dirac field can constitute dark matter. It is stable and electrically neutral. Since the 
model we consider has two components of dark matter, the total relic density Vth 2 will be shared by 
the neutralino and Dirac particles. In the analysis we assume that the dark matter densities g D , g x 
for the two components in the galaxy are proportional to their respective relic densities such that 
sum is the total cold dark matter (CDM) density 



— n CDM h 2 = n x h 2 +n D h 2 . (59) 

Q\ X (Majorana) (Dirac) 



The annihilation cross section of DD into quarks and leptons via the Z', Z" poles is given by 

&DD->ff = Adj\Pz> — Pz>>\ 2 , (60) 
where the poles and couplings enter as 

P v = {s-M 2 +iT v M v )-\ V = (Z',Z"), (61) 

A °-< = itr<< 2M ° + s)i2m > + s) V <62) 

9dj = 9blQblj9xQx sin 29bl, (63) 
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and where s = 4M£/(1 - v 2 /4), 6 = Q(s - Amj), and AT/ = (1, 3) for (leptons, quarks). The 
relevant partial Z', Z" decay widths were given in Table (1) . In addition the Z', Z" can decay into 
the Dirac sector: 

r R M z ,gl( 2Ml\( AMlV' 2 

r ™ = -^rl v 1 + -mi) y- mi) ' (64) 

where = Q(Mz> — 2M D ) and g D = gxQx cos9 BL . The partial decay width of the Z" is ob- 
tained with M Z t — > Mz" and cos9 B l — > sin-9 BL in Eq.(l64i The relic density can be calculated 



990. The relic 



by integration over the poles. For the technique of integrating over a pole see 1197 - 

density for the 2 components of dark matter can be calculated 07I1 where for the Dirac component 

r x F,D ^ 1 07 x 10 9 GeV -1 

Vt D h 2 = C D J D \ J D = V {crv} DD ^ ff dx, C D = 2x — — . (65) 

Jo j v9 -^pi 

In Fig. © we exhibit a satisfaction of the relic density within the WMAP constraint so that 

^Dirac = ^ * 1/2, (66) 

where the black bands in Fig. (J6]) show a presumed fraction of the the total relic abundance. 

Now, unlike the cases studied previously with the Stueckelberg mass growth, here the dark 
Dirac fermion does not carry a milli-charge and is electrically neutral. The reason the Dirac 
fermion is neutral is because there is no mixing of the Stueckelberg gauge field with the hyper- 
charge vector boson. Because of its electrical neutrality and unlike a milli charged particle it cannot 
be stopped by the atmosphere or by dirt and rock in the Earth before it reaches an underground 
detector. The effective Lagrangian describing the scattering of a dark Dirac fermion from a quark, 
in the limit of low momentum transfer, is given by £ eff = C^D^Dq^^q. The corresponding spin 
independent D-proton cross section is 

SI _ ^Dp r 2 ( 1 j 1 2 A (f<1 ^ 

a ° p ~ vr \M Z , + Mf„ Mf,M|J ' ^ ] 

where G = g B L sin BL g x Qx cos 9 BL and fx Dp is the reduced mass. 

Interestingly, for mixing of the size considered in Fig. (Q]), (sin9 B L £ [0.01,0.05]) and for 

natural size couplings gx = gBL = O(gy) and Qx = ±1 one obtains a spin independent cross 

sections which are of the size 

erg, ~ 10- 45±1 cm 2 , M z < ~ (200 - 300) GeV. (68) 

26 



Z prime Decays into SM Z prime Decays into HS 




FIG. 6: An exhibition of the relic density of the Dirac component of dark matter for various values of •flfj} rac 
which is the ratio the Dirac dark matter mass to the Stueckelberg Z' mass. The black bands represent about 
half the relic abundance. For the analysis we fix gBL = 0.35, gx = 0.1, Qx = 0.5 . The (blue/darker) 
curves have the Z' mass running in the range 200-500 GeV in steps of 100 GeV. We note that for fixed 
couplings, as Mz> gets heavier the curves become more narrow. The (magenta/lighter) curves correspond 
to Mz> = 250 GeV with 9bl = (0.02 — 0.05). Similarly, as 9bl becomes progressively smaller for 
otherwise fixed couplings and fixed Z' mass, the curves become more narrow. The right panel is the case 
when the Z 1 decays mostly into the hidden sector Dirac fermions, i.e., it is the case where Z 1 — > DD is 
kinematically allowed and in this case the dileptonic signals at the LHC will be depleted. The left panel is 
the case where Z' — > DD is kinematically disallowed and in this case the Z' will decay exclusively into the 
SM particles and thus the dileptonic signal from the process pp — > Z' — > l + l~ will be visible. 



Since M D ^> m p , fx Dp ~ m p , af^ p is essentially independent of M D . However, compatibility 
with the WMAP data for the thermal relic density, restricts the ratio R^ irac — 1/2. Using this 
constraint the spin independent cross section for the case M|„ 3> M§, is given by 



y 71 M|, l07T 

which now has a very strong dependence on the Dirac mass. The numerical size of as a 
function of the Dirac mass is exhibited in Fig. (5), and the analysis shows that the afj p predicted 
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FIG. 7: Illustrative curves. At any given point on this plot there exists a funnel where the relic density can 
be satisfied for perturbative size coupling via the relic density invariant -R^irac- The particular values of the 
parameters on thes curves are (9bl = -001, gx = 0.1,Qx = 1/2), {@bl = 0.01,gx = 0.5, Qx = 1), and 
(9 BL = 0.05, g x = 1/2, Qx = 1), where g BL = 0.35 and R%\ iac ~ 1/2. 

by the model is accessible in the XENON experiment. In fact for given values of Qbl-, &bl, QxQx 
the current limits from XENON 100 already put lower limits on the Dirac mass. We can also use 
the current upper limit on a SI from the XENON100 experiment which gives o SI = 7 x 10~ 45 cm 2 
for a WIMP mass of 50 GeV, to put a general constraint on \G\/M D . so that 

\G\/M 2 D < 3 x 10~ 8 (M D in GeV). (70) 

We note again that the preceding analysis is very different from the previous Stueckelberg analyses 
where the Dirac fermion in the hidden sector develops a milli charge. As already pointed out this 
arises in models where one mixes the Stueckelberg gauge boson with the hypercharge gauge field. 
In this case the scattering of the Dirac fermion from a quark will have not only the Z' pole in the 
t-channel but also a Z boson pole and a photon pole as well. In the present model the Z and the 
photon pole are both absent. The Dirac dark matter candidate is electrically neutral. 

As mentioned earlier, for Mz> ~ 2Md, the relic density will always be satisfied for perturbative 
size couplings. For M Z ' < 2M D but close to 2M D the Z' signal will manifest at colliders and the 
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M z > <% {Obl = (0.03)) a% p (6 BL = (0.06)) 



GeV 



cm 



cm 



200 1.9xl0~ 44 

300 3.7 xlO" 45 

400 1.2xl0~ 45 

500 4.8 xlO" 46 



7.5 xlO- 44 
1.5 xlO- 44 
4.7 xlO" 45 
1.9 xlO" 45 



TABLE II: Approximate values of the spin independent scattering cross section for the Dirac component 
of dark matter for sample models. The second and third columns have Obl = (0.03,0.06) respectively. 
The first row is on the edge of the discovery limits from the both XENON and the Tevatron data and is 
being probed by the LHC. For a given dark matter mass Mz> ~ 2Mjj in order to satisfy (f2/i 2 )\yMAP- 
Model parameters are otherwise fixed as in Figure(4). The middle column of this table corresponds to the 
blue/dark curves in Fig. (7), while the magenta/light region is found to be constrained by the XENON data. 
Models consistent with the relic density constraint and the XENON constraint are therefore favored if the 
relic density is satisfied closer to the pole which is obtained for relatively smaller coupling and/or larger 
M Z : 



relic density can also be satisfied. However, for the case M Z ' > 2Mp, while the relic density 
can be satisfied, the Z' signal becomes suppressed due to the branching ratio into the hidden 
sector o verta king the branching ratio in the visible sector in the presence of mass and kinetic 
mixings 1 10011 . In addition, the Breit-Wigner enhancement of the annihilation of Dirac particles in 
the halo [41] can be operative very close to the pole and the following three possibilities become 
simultaneous observables: 

1 . Observation of a very light and narrow Z' vector boson in the dilepton channel at the LHC 
(see also |5|]). 



2. Observation of the flux of positrons via Satellite data (PAMELA/FERMI) 110 



Tom the 



Breit-Wigner Enhancement in the dark matter annihilations in the galactic halo 114 ill consis- 
tent with WMAP data D49Q. 



3. Relic abundance of dark matter split between a neutralino and dark Dirac (see also 1371]) 
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4. Observational prospects for the corresponding Dark Dirac component in direct detections 
experiments such as XENON (analyzed here for the neutral dark Dirac particle via the 
Stueckelberg mechanism). 

Let us add , that just recently, the 730 kg days of the CRESST-II Dark Matter Search was 
released fl 102TI . Two preferred regions are reported on, and one such region appears close to the 
CoGeNT preferred region [11030. Very low mass neutralino dark matter with MSSM field content 
and cross s ection s of the size needed to explain the CoGeNT are not consistent with the collider 
constraints [104]. T his r esult has been confirmed by the LHC with its updated constraints on the 
SUSY Higgs sector Iil05ll . wherein large tan (3 and low mass SUSY Higgs of the size needed to 
explain the spin-independant scattering are further excluded. The preferred region reported by 
CRESST-II with heavier dark matter mass may be accommodated for a thermal relic with relic 
density satisfied via the Z-pole in the MSSM. Such could arise with non-universal gaugino masses 
at the the high-scale (see [93]) leading to WIMP masses close to 45 GeV. The far boundary of the 
CRESST-II 2er region terminating clos e to 5 5 GeV may also be achieved with relic density satisfied 



via the Higgs pole (see the analysis of [ 



SUSY Higgs sector from the LHC 



106D). A dedicated analyses with the new constraints on the 



1051 would be needed to make a more definitive statement 



however the CRESST-II results at these potential dark mater masses do not correspond to reported 



event rates with CDMS or XENON [50, 



5lh . The extended model class we discuss can produce 



spin independent cross sections with larger cross sections than that of the MSSM via the Dirac 
component of Dark Matter (see Fig.©). 



VII. DISCRIMINATING STUECKELBERG FROM MODELS WITH SPONTANEOUS BREAK- 
ING 

One may discriminate between the Stueckelberg mass growth for a B — L gauge boson in the 
models discussed here and other models where the mass growth for the B — L gauge boson occurs 
by spontaneous breaking. In the above, we have already discussed the mass growth of a B — L 
gauge boson by the Stueckelberg mechanism. For the case when the mass growth occurs via 
spontaneous breaking there are two possibilities: (i) spontaneous symmetry breaking of U(\)b-l 
occurs violating R-parity invariance, (ii) spontaneous symmetry breaking of U(1)b-l occurs 
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without violating R-parity invariance. We discuss these two cases below individually. 



A. Spontaneous Symmetry Breaking of B — L and R-parity Violation 

The simplest example of this is when we consider the superpotential of Eq.(fl"3l). Let us 
assume that the potential of the v c field is such that it develops a VEV. In this case one will 
have a spontaneous breaking of not only B — L but also of R-parity as indicated by the term 
LH u (u c ) in Eq.(fT3~l) after v c develops a VEV. In the mass diagonal basis it will lead to other 
R-parity violating terms, i.e., LLe c and QLd c . Here the LSP is no longer stable and specifically 
the neutralino cannot be a dark matter particle. Further, since the neutralino is not stable, the 
signals of supersymmetry for this case will be very different at hadron colliders. Specifically if 
the neutralino decays inside the detector, there will be no missing energy signatures which are 
the typical hallmarks of supersymmetry signatures with R-parity symmetry. Further, for the case 
when there is a spontaneous breaking of R-parity symmetry via the VEV growth of the right 
handed sneutrino, there will be D term contributions to the slepton squared masses proportional 
to g BL {v c ) 2 \m. Such terms are absent for the case when the mass growth for the B — L gauge 
boson occurs preserving R-parity invariance as discussed below. 



B. B — L Models for R-parity Conservation 

We further consider now the possibility that B — L symmetry is broken but a residual R- parity 



symmetr y stil l persists. This is indeed possible following the general line of reasoning of F107I1 



(see also H108I1 '). Thus consider additional fields in the theory such as a vector like multiplet which 



has the SU (3)e ® SU(2) L £g> U (l)y £g> U (1) b-l quantum numbers as follows 

$ ~ (1, 1, 0, Qbl) i $ ~ (1, 1, 0, Q BL ). (71) 

Let us suppose that one manufactures a potential so that VEV formation for the fields $ and $ 
occurs. In this case B — L will be broken. However, as long as 3(5 — L) is an even integer R-parity 
will be preserved. This means that the residual theory will have a Z 2 R-parity symmetry. Thus, for 
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example, the VEV formation of a scalar field with 3(B — L) = ±2 will violate B — L but preserve 
R-parity. In the process of the mass mass growth of the B — L gauge boson, one combination of 
the imaginary parts of $° and $° will be absorbed while there would three spin zero fields: 2 CP 
even and one CP odd (the part orthogonal to the imaginary parts of $° and $° which is absorbed) 
Higgs field. In contrast for the £7(1) b-l ® U(l)x model discussed here, one is left with only two 
additional scalars, px,Pbl, or p,p', which are both CP even. Specifically there is no additional 
CP odd Higgs boson for the Stueckelberg models. So this provides a discrimination between the 
two models. 

There are several interesting and distinguishing features between the U(l) b—l <E> U{l) x model 
and the U(\)b-l model. This difference can be seen by comparing Eq.(|23l vs Eq.(|30l). Thus in 
Eq.(l23l one finds that the mass growth ofaB — L gauge boson by spontaneous breaking or by the 
Stueckelberg mechanism would require the gauge boson to be very heavy. Thus for gnj, ~ 1, one 



will typically have a mass of the B — L gauge boson to be greater than ~ 6 TeV [69, 



109]. In 



contrast, from Eqs.(l3~0l) we find that in the U(l)x <8> U(l) b-l model, there are two extra massive 
gauge bosons beyond what one has in the Standard Model. Thus the heavier one, i.e., the Z" gauge 
boson, is indeed several TeV in mass. However, the Z' boson we discuss can be much lighter, and 
can lie in the few hundred GeV range. Thus the observation of a low lying Z' with decay branching 
ratios characteristic of a B — L gauge boson will be a clear indication of the Stueckelberg model 
involving mixing of U(l)x and U (I)b-l discussed here. 



VIII. CONCLUSION 



In this work we have proposed the Stueckelberg mechanism for the mass growth of a B — L 
gauge boson. It was then shown that under the constraints of charge conservation and the absence 
of a Fayet-Iliopoulos D term, that R-parity cannot be spontaneously broken in the minimal model 
of radiative electroweak symmetry breaking. The above is in contrast to models where the mass of 
the B — L gauge boson is generated by the Higgs mechanism through the VEV formation for the 
field v c which breaks R-parity. 

A comparison to the case where the B — L symmetry is spontaneously broken but the R-parity 
symmetry is preserved was also given and its distinguishing features from the Stueckelberg mass 
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growth for the B — L gauge boson are uncovered. Further, we analyzed a U(l)x <8> U(1)b-l 
Stueckelberg extension of MSSM where a massive Z' boson with B — L interactions can lie in 
the sub TeV region, i.e, M z > < 1 TeV. The observation of a Z' in the sub TeV region with B — L 
quantum numbers deduced via branching ratios into charged leptons will provide a test of the 
U(l)x ® U(l) b-l Stueckelberg extension discussed here. 

Other tests of the proposed Stueckelberg models were also discussed. This includes an analysis 
of the production and decay of the Stueckelberg spin boson p which has only loop decays into 
SM final states via sfermion loops. An interesting decay of the p is into 77 was analyzed and 
shown to have the possibility of observation at the LHC with y/s = 14 TeV. 

With hidden sector Dirac fermions in the U(l)x <S> U(\)b-l Stueckelberg extension, two 
component dark matter manifests, with one component being either the MSSM neutralino or the 
Stueckelberg neutralino and the other component being a neutral Dirac fermion. An analysis 
of the relic density for the Stueckelberg neutralino and the Stueckelberg neutralino-proton spin 
independent cross section were also discussed. An analysis of the second dark matter component 
consisting of the Dirac fermion as dark matter was also given and it was shown that the current 
XENON100 data already puts constraints on the Dirac fermion mass and mixing angles. The 
constraints from the XENON 100 data and the LHC data on the couplings of the Z 1 boson and 
dark Dirac fermion were shown to be comparable, both of which limit the mixing of the B — L 
and dark sector. Thus the proposed model produces LHC and dark matter signals at mass scales 
that are accessible to such experiments and will be tested further as the new data comes in. 
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